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(54) Method and apparatus for reducing three-dimensional shape data 



(57) A method and apparatus for reducing three-di- 
mensional shape data, includes calculating estimation 
values for portions to be deformed by converging two or 
more vertices (Vj) of a polygon model (M|), and reducing 
the number of data for the polygon model (Mj) by con- 
verging two or more vertices (Vj) of the polygon model 
(Mj) when the estimation values are equal to or below 
the predetermined permissible value. Estimation values 



for surfaces to be deformed by shrinking edges or sur- 
faces are calculated based on distances between the 
respective surfaces and all the original vertices (Vj) in- 
volved in the surface deformation. Before each data re- 
duction, the portion having been involved in all the pre- 
vious data reductions is defined as a reduction prohibi- 
tion area, and a succeeding data reduction is applied to 
a portion other than the reduction prohibition area. 
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Description 

[0001] This invention relates to a metliod and apparatus for reducing tliree-dimensional sliape data forming a polygon 
model. 

5 [0002] A huge number of measurement data are obtained by, for example, a three-dimensional shape measuring 
apparatus having quite a high space resolution if an image quality is high. Accordingly, if the data have a higher image 
quality than necessary, it disadvantageously results in a large burden on a storage device during succeeding processing 
and a considerable reduction in processing speed. Thus, it is necessary to ease this problem as much as possible by 
reducing the number of measurement data. However, in such a case, if a difference between a polygon model (here- 
to inafter, referred to as approximated polygon model) formed by the measurement data after the data reduction and a 
polygon model formed by the actually measured three-dimensional shape data (hereinafter, referred to as original 
polygon model) is excessively large, the measurement data after the data reduction cannot be effectively utilized. 
Therefore, it is necessary to control error between the approximated polygon model and the original polygon model in 
the data reduction. 

15 [0003] Particularly in the case of three-dimensional shape data obtained by a three-dimensional shape measuring 
apparatus, actually measured three-dimensional shape data are converted into an approximate value by the data 
reduction. Thus, error between the approximated polygon model and the original polygon model is preferably controlled 
using the actually measured three-dimensional shape data as a reference. 

[0004] On the other hand, various methods have been proposed to reduce the number of three-dimensional data 
20 forming a triangular polygon model by approximating a model within a desired permissible error range. 

[0005] For example, In "Hierarchical Approximation of Three-Dimensional Polygon IVIodels" (by Horil^awa and Tot- 
suka, pp. 3-7 of Proceedings of the 5^^ Sony Research Forum, 1995) is disclosed a method for reducing the number 
of three-dimensional shape data by estimating contributions of respective edges forming a triangular polygon model 
in the form of an estimation function expressing a volume which changes in the case that edges are deleted, and 
25 deleting edges, for example, in the case that they are parallel to surrounding surfaces, surfaces are small and/or edges 
are short. 

[0006] Further, in "Simplification Envelopes" (by Cohen Jonathan, Amitabh Varshney, Dinesh Manocha, Greg Turk, 
Hans Weber, Pankai Agarwal, Frederick Brooks and William Wright, pp. 119-128 of Proceedings of SIGGRAPH 96 in 
Computer Graphics Proceedings, Annual Conference Series, 1 996, ACM SIGGRAPH) is disclosed a method for gen- 
30 erating models obtained by reducing an original polygon model by an error e and enlarging it by the error e and reducing 
the number of data of the original polygon model such that an approximate polygon model lies between the enlarged 
and reduced models. 

[0007] Further, in "Surface Simplification Using Quadric Error Metrics" (by Michael Garland and Paul S. Heckbert, 
Carnegie Mellon University, Proceedings of SIGGRAPH 97, 1997) is disclosed a method for reducing the number of 

35 three-dimensional shape data forming a triangular polygon model by converging two or more vertices (three-dimen- 
sional shape data) forming an edge or surface of the triangular polygon model to one vertex. In this method, a distance 
between the converged vertex due to edge shrinking or surface shrinking and the respective planes forming the original 
triangular polygon model influenced by the edge shrinking or surface shrinking is controlled as an error in the data 
reduction. The number of data is reduced within a specified permissible error range by performing such an edge shrink- 

40 ing or surface shrinking that this error lies within a specified permissible range. 

[0008] In this method, specifically, when the edge shrinking is applied for one edge, the error estimation is executed 
for the other edges or surfaces of the polygon model that are to be influenced by the edge shrinking because the 
shrinking of the one edge changes edges and surfaces surrounding the one edge. After the error calculation is com- 
pleted, the shrinking of the one edge is determined. This processing is repeated to reduce the number of the three- 

45 dimensional shape data within a predetermined permissible range. 

[0009] Since the data reducing method disclosed in "Hierarchical Approximation of Three-Dimensional Polygon Mod- 
els" only deletes edges based on evaluation values thereof, an error in the approximation of a model by the edge 
deletion cannot be estimated although the three-dimensional data of characteristic portions of a polygon model remain. 
Thus, this method is effective in the case of reducing the number of data of a polygon model for merely displaying a 

50 three-dimensional image such as virtual reality and games in computer graphics. However, this method cannot be 
effective in the case that the error between the original measurement data and the data of the approximated model 
after the data reduction needs to be controlled. 

[0010] On the other hand, the data reducing method disclosed in "Simplification Envelopes" can strictly control such 
error. However, it is necessary to generate polygon models by enlarging and reducing the original polygon model by 
55 an error e and save the data of these polygon models during the data reduction. Thus, the number of data relating to 
the processing becomes huge, causing the problem of large burden on the memory capacity and calculations for the 
data reduction. 

[0011] Further, the method disclosed in "Surface Simplification Using Quadric Error Metrics" controls the distances 
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from the vertices of the polygon model after the approximation from the surfaces forming the original polygon model 
as errors. Thus, this method is not necessarily a suitable method in the case that data reduction is desired to be 
performed by controlling errors of the polygon models after the approximation from the vertices forming the original 
polygon model like the measurement data of the three-dimensional shape measuring apparatus. 
5 [0012] Further, in the case of measurement data obtained by measuring rough surfaces or cut surfaces of a meas- 
urement object by means of a three-dimensional shape measuring apparatus, e.g., in the case that the orientation of 
surfaces locally suddenly change in a zigzag manner in a polygon model, errors for these zigzag surfaces are estimated. 
Therefore, errors are estimated as large values as the data reduction progresses, thereby presenting the problem of 
insufficient data reduction. 

10 [0013] Furthermore, each time the edge or surface shrinking is performed, it is necessary to estimate errors with 

respect to the edges or surfaces of the shrunk polygon model that are to be influenced by the edge or surface shrinking. 

This will increase the number of calculations, consequently causing the problem of hindering the data reduction. 

[0014] It is an object of the present invention to provide a method and apparatus for reducing three-dimensional 

shape data which are free from the problems residing in the prior art. 
15 [0015] It is another object of the present invention to provide a three-dimensional shape data reducing method and 

apparatus which enable higher speed processing. 

[0016] It is another object of the present invention is to provide a three-dimensional shape data reducing method 
and apparatus which enable proper data reduction based on precise error estimation for measurement data obtained 
by a three-dimensional shape measuring apparatus. 
20 [0017] According to an aspect of the invention, a method for reducing three-dimensional shape data, comprises 
calculating estimation values for portions to be deformed by converging two or more vertices of a polygon model, and 
reducing the number of data for the polygon model by converging two or more vertices of the polygon model when the 
estimation values are equal to or below the predetermined permissible value. 

[0018] The portions to be deformed may be surfaces of the polygon model. Estimation values for surfaces to be 
25 deformed by shrinking edges or surfaces may be calculated based on distances between the respective surfaces and 

all the original vertices involved in the surface deformation. 

[0019] Alternatively, the number of data for the polygon model may be reduced by converging two or more vertices 
of one portion of the polygon model after another portion, repeatedly. Before each data reduction, the portion having 
been involved in all the previous data reductions is defined as a reduction prohibition area, and a succeeding data 

30 reduction is applied to a portion other than the reduction prohibition area. 

[0020] According to another aspect of the invention, an apparatus for reducing three-dimensional shape data, com- 
prises an estimation value calculator for calculating estimation values for portions to be deformed by converging two 
or more vertices of a polygon model, and a data reducing device for reducing the number of data for the polygon model 
by converging two or more vertices of the polygon model when the estimation values are equal to or below the prede- 

35 termined permissible value. 

[0021] The estimation value calculator may calculate estimation values for surfaces to be deformed by shrinking 
edges or surfaces based on distances between the respective surfaces and all the original vertices involved in the 
surface deformation. 

[0022] Also, the data reducing device may reduce the number of data for the polygon model by converging two or 
40 more vertices of one portion of the polygon model after another portion, repeatedly. Further, there may be provided a 
judging device forjudging based on a calculated estimated value whether a portion of the polygon model is permissible 
for the deformation, and a prohibition area defining device for defining a portion having been involved in all the previous 
data reductions executed by the data reducing device as a prohibition area to keep a succeeding data reduction from 
being applied. 

45 [0023] These and other objects, features and advantages of the present invention will become more apparent upon 
a reading of the following detailed description and accompanying drawings, in which: 

FIG. 1 is a block diagram showing a construction of a three-dimensional shape measurement data processing 
system according to an embodiment of the invention; 
50 FIG. 2 is a block diagram showing a construction of a shape data processor; 

FIG. 3 is a diagram showing an example of approximation of a polygon model by data reduction; 
FIG. 4 is a concept diagram showing a data reducing method by edge shrinking; 
FIG. 5 is a concept diagram showing a data reducing method by surface shrinking; 

FIG. 6 is a diagram showing an estimation value used in a data reducing method according to the embodiment; 
55 FIG. 7 is a flowchart showing a procedure of a three-dimensional shape data reduction processing; 

FIG. 8 is a diagram showing edge shrinking in a portion where a discontinuous surface is formed on a polygon 
model by a random noise data; 

FIG. 9 is a diagram showing a difference between an error estimation value according to the embodiment and a 
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conventional error estimation value in the data reduction shown in FIG. 8; 

FIG. 1 0 is a diagram showing shrinking of a zigzag surface of a polygon model having such a surface; 
FIG. 11 is a diagram showing a conventional error estimation value in the data reduction shown in FIG. 10; 
FIG. 12 is a diagram showing an error estimation value according to the embodiment in the data reduction shown 
5 in FIG. 10; 

FIG. 13 is a flowchart showing another embodiment of the three-dimensional shape data reduction procedure; 
FIG. 14 is a diagram showing an example of a data reducing process by the operation procedure shown in FIG. 
13; and 

FIG. 15 is a diagram showing a modification of the flowchart shown in FIG. 13. 

10 

[0024] Referring to FIG. 1, a three-dimensional shape measurement data processing system (hereinafter, "shape 
data processing system") is provided with a three-dimensional shape measuring apparatus 1 for obtaining three-di- 
mensional shape data of a measurement object G to be measured, an external storage device 2 for storing three- 
dimensional shape data obtained by the three-dimensional shape measuring apparatus 1 , a shape data processor 3 

15 for reducing the three-dimensional shape data obtained by the three-dimensional shape measuring apparatus 1 by a 
data reducing method to perform a specified data processing, an external storage device 4 for storing the three-di- 
mensional shape data processed by the shape data processor 3, a display device 5 for displaying a three-dimensional 
image of the measurement object G such as a polygon model using the three-dimensional shape data after data 
reduction, and an external processing device 6 such as a CAD/CAM (computer aided design and computer aided 

20 manufacturing) device in which the three-dimensional shape data processed by the shape data processor 3 are utilized. 
[0025] The three-dimensional shape measuring apparatus 1 measures a surface shape of the measurement object 
G using a conjugated autofocus system. As shown in FIG. 2, the three-dimensional shape measuring apparatus 1 
basically includes an illuminating unit 1 01 for generating illumination light, a confocal optical system 1 02 for projecting 
illumination light from the illuminating unit 101 onto the measurement object G and introducing the illumination light 

25 reflected by the measurement object G to a photodetector 1 03, and the photodetector 1 03 for receiving the illumination 
light reflected by the measurement object G and outputting it after photoelectrically converting it into an electrical signal, 
and a beam splitter 104 for splitting the light reflected by the measurement object G and having transmitted through 
the confocal optical system 102 (hereinafter, referred to as merely "optical system 102"), and a calculator 105 for 
calculating coordinates of surface positions of the measurement object G using light reception signals outputted from 

30 the photodetector 1 03. 

[0026] The three-dimensional shape of the measurement object G is measured by placing the measurement object 
G on a measurement table 1 06 provided in front of the three-dimensional shape measuring apparatus 1 and measuring 
the surface shape of the measurement object G at each height position while elevating and lowering the measurement 
table 106 at specified intervals along height direction. Specifically, assuming that Y-, Z- and X-directions denote the 

35 height direction of the measurement object G, the optical axis direction of the optical system 1 02 and a direction normal 
to Y-direction and Z-direction, respectively, the three-dimensional shape of the measurement object G is measured by 
measuring Z-coordinate of the measurement object G at each Y-coordinate while changing the Y-coordinate of the 
measurement object G at specified intervals. The three-dimensional shape data (x, y, z) of the measurement object G 
measured by the three-dimensional shape measuring apparatus 1 are inputted to the shape data processor 3. 

40 [0027] The three-dimensional shape measuring apparatus 1 is provided with a photoelectrically converting device 
1 03a arranged in parallel to X-axis and formed of a line sensor, and the coordinate (Z-coordinate) of a surface position 
of the measurement object G is measured using a light reception signal for each pixel (each pixel position corresponds 
to an X-coordinate) of the photoelectrically converting device 103a. An illumination light in the form of a slit extending 
in X-direction is emitted from the illuminating unit 1 01 so that the light reflected by the measurement object G is incident 

45 on the entire photoelectrically converting device 103a. 

[0028] The Z-coordinates of the measurement object G are measured as follows. First, illumination light in the form 
of the slit extending in X-direction is emitted from the illuminating unit 101 provided at an image-side focusing position 
of the optical system 102. This illumination light is gathered at an object-side focusing position of the optical system 
1 02 to form a line beam extending in X-direction at this focusing position. Subsequently, the photoelectrically converting 

50 device 103a is repeatedly exposed in a specific cycle while the object-side focusing position (i.e., line beam) of the 
optical system 102 is moved in the optical axis direction until this focusing position reaches the inside of the measure- 
ment object G. An exposure position where the level of the light reception signal is at maximum is calculated by the 
calculator 105 for each pixel of the photoelectrically converting device 103a using the levels of the light reception 
signals obtained at the each exposure position, and this position is assumed to be the Z-coordinate of the outer surface 

55 of the measurement object G. 

[0029] The shape data processor 3 is provided with a polygon generator 31 , a data reducing device 32, a permissible 
value input device 33 and a data output device 34. The polygon generator 31 generates a polygon model of the meas- 
urement object G based on the three-dimensional shape data (x, y, z) of the measurement object G inputted from the 
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three-dimensional sliape measuring apparatus 1. The data reducing device 32 reduces the three-dimensional shape 
data by edge shrinking or surface shrinking of the polygon model. In other words, as shown in FIG. 3, a precise polygon 
model A based on the measurement data is approximated to a polygon model B within a specified permissible error 
range by suitably reducing the vertices (three-dimensional shape data) of the polygon. 

[0030] The data reducing device 32 includes an estimation value calculator 321 , a estimation value sorter 322 and 
an edge/surface shrinking section 323 to carry out a data reduction processing to be described later. 
[0031 ] The estimation value calculator 321 calculates estimation values used to judge whether the respective edges 
or surfaces can be shrunk before shrinking edges or surfaces forming the polygon model by converging two or more 
vertices to one vertex. This estimation value is described later. 

[0032] The estimation value sorter 322 includes a calculator for calculating a minimum estimation value by rearrang- 
ing a plurality of estimation values calculated for the respective edges or surfaces by the estimation value calculator 
321 in an increasing or decreasing order. The edge/surface shrinking section 323 compares the minimum estimation 
value calculated by the calculator of the estimation value sorter 322 with a permissible value inputted via the permissible 
value input device 33 and performs an edge or surface shrinking corresponding to the minimum estimation value when 
the minimum estimation value is equal to or below the permissible value. The estimation values and an edge shrinking 
method are described later. 

[0033] The permissible value input device 33 is operated by an operator to input the permissible value used to judge 
whether the edges or surfaces can be shrunk as described above. This permissible value corresponds to a permissible 
range within which the edge or surface shrinking is assumed to have a little influence on characterizing portion(s) of 
the polygon model in approximating the polygon model, and is suitably inputted by the operator according to the shape 
of the measurement object G. 

[0034] The data output device 34 is an interface for outputting the reduced three-dimensional shape data (three- 
dimensional shape data forming the approximated polygon data) to the external storage device 4, the display device 
5 and the external processing device 6. 

[0035] At least part of devices of the shape data processor 3 may be provided in another independent processing 

unit, and their functions may be realized by cooperation of the processing units and the shape data processor 3. 
Alternatively, these devices may be processing units which are separately provided in the shape data processor 3 and 
cooperate with each other. Conversely, at least two functions of these devices may be realized by a processing program 
of one processor, e.g., the shape data processor 3 or a processor provided therein. 

[0036] The external storage device 4 is adapted to store the data in an external storage medium such as a magnetic 
disk, an optical disk or a magnetooptical disk. Upon input of the three-dimensional shape data from the shape data 
processor 3, the external storage device 4 stores the data in an external storage medium in a specified file format. 
[0037] The display device 5 is an electronic display device such as a cathode-ray tube or a liquid crystal display 
device. Upon input of the three-dimensional shape data from the shape data processor 3, an image of the approximated 
polygon model of the measurement object G is displayed on the display device 5 based on the inputted data (see the 
polygon model B of FIG. 3). 

[0038] The external processing device 6 is a device for fabricating, for example, a mold or a model having substantially 
the same configuration as the surface configuration of the measurement object G. Upon input of the three-dimensional 
shape data from the shape data processor 3, a mold or the like is automatically fabricated by cutting a work based on 
the inputted three-dimensional shape data. 

[0039] Referring back to FIG. 1 , at least part of the measuring apparatus 1 and the devices 2, 4, 5, 6 may be incor- 
porated in the shape data processor 3. 

[0040] Referring to FIGS. 4 and 5, next, the data reducing method will be described. The data reducing method by 

edge shrinking is such a method for reducing the number of data that, in the case that a polygon model is defined by 
vertices v^ to v^q shown in FIG. 4, two vertices v^, are, for example, so converged to a vertex v-|' as to shrink an 
edge Ed-, 2 connecting the vertices v-,, V2. 

[0041] In the data reducing method by edge shrinking (hereinafter, referred to as "edge shrinking method"), estimation 
values for estimating the influence of the edge shrinking on a change of the configuration of the polygon model for all 
edges Edy (i, j = 1, 2, ...10, combinations of vertices adjacent to each other when i^^j) are first calculated, and the 
vertices Vj, Vj at the opposite ends of the edge Edy corresponding to a minimum estimation value are converged to a 
specified vertex Vj' for data reduction if the minimum estimation value is equal to or below a permissible value. 
[0042] Subsequently, the edge shrinking is performed in a similar manner for the polygon model (approximated 
polygon model) in which the vertices Vj, Vj are converted to the vertex Vj', thereby further reducing the data. Thereafter, 
the edge shrinking is performed edge by edge in a similar method until no more edge shrinking can be performed for 
the polygon model after edge shrinking (i.e., the approximation of the polygon model is repeated until the minimum 
estimation value becomes larger than the permissible value) to reduce the data. 

[0043] It should be noted that the vertex vf is an approximated point set in a specified position between the vertices 
Vj, Vj such as a middle point of the edge Edy as disclosed in the aforementioned paper "Hierarchical Approximation of 
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Three-Dimensional Polygon Models". 

[0044] The data reducing method by surface shrinking (hereinafter, referred to as merely "surface shrinking method") 
is such a method for reducing the number of data that, in the case that a polygon model is defined by vertices to 
v-11 as shown in FIG. 5, vertices v^ , V2, V3 are so converged to a vertex v^' as to delete a surface P^23 enclosed by the 
vertices v-,, V2, V3. 

[0045] Basically, the data reducing method by surface shrinking is carried out in a procedure similar to that of the 
data reducing method by edge shrinking. Specifically, estimation values for estimating the influence of the surface 
shrinking on a change of the configuration of the polygon model for all surfaces Py^ (i, j, k = 1 , 2, ...10, combinations 
of vertices adjacent to each other when i^j^^^k) are first calculated, and the vertices Vj, Vj, v^^ defining the surface P^^^ 
corresponding to a minimum estimation value are converged to a specified vertex Vj' for data reduction if the minimum 
estimation value is equal to or below a permissible value. Thereafter, the surface shrinking is performed surface by 
surface in a similar method until no more surface shrinking can be performed for the polygon model after surface 
shrinking to reduce the data. 

[0046] In the case that the three-dimensional shape data are successively reduced by the edge reducing method or 
surface reducing method, the three-dimensional shape data forming the measurement object G after the data reduction 
processing are different from those of the measurement object G actually measured by the three-dimensional shape 
measuring apparatus 1 and serve as approximated three-dimensional shape data. Thus, it is necessary to perform the 
data reduction processing while controlling an error of the polygon model after the data reduction (hereinafter, referred 
to as "approximated polygon model") from the polygon model formed by the actually measured three-dimensional 
shape data (hereinafter, referred to as "original polygon model"). 

[0047] In the data reducing method according to this embodiment, if it is assumed that a triangular polygon P^b^ is 
approximated to a triangular polygon P^^^^' by, for example, approximating vertices a, b, c (actually measured three- 
dimensional shape data) to vertices a', b', c' as shown in FIG. 6 through, for example, repetition of the edge shrinking, 
distances of a surface including the triangular polygon P^^^' after the approximation from the original vertices a, b, c 
involved in the approximation of this surface are controlled as errors. 

[0048] In the example of FIG. 6, the vertices involved in the approximation are assumed to be the vertices a, b, c in 
order to simplify the description. However, in an actual processing, vertices other than the vertex a are involved in the 
approximation, for example, while the vertex a is converged to the vertex a'. Therefore, the above errors include dis- 
tances of these vertices from the triangular polygon Pabc - 

[0049] Next, the estimation value of the error used in the data reducing method according to this embodiment is 
described. Generally, a surface p in a three-dimensional space is expressed by Equation (1): 



[0050] On the other hand, if (v^, vy, v^) denotes the coordinates of point v in the three-dimensional space, a distance 
d between the point v and a surface p is calculated by Equation (2): 



A- x+ B- y+ C- z+ D = 0 



(1) 



where 



cl = A -v^ + B -v^ + C -v^ + D 



(2) 



and can also be expressed by matrices as in Equation (3): 



V 



(3) 



where 
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a 
b 
c 
d 



P^is a transposed matrix of P. 

[0051] Accordingly, the square of the distance d is expressed by Equation (4): 

af^ = (P^- vf = (P^- V)-(V^ • P) = -{V- V^) - P= P^ • Kv ■ P 
where is a transposed matrix of V 



(4) 



V V V V 



V V V 

^ y y 

V V V V 



V V 



[0052] Equation (4) expresses the square of the distance d between one plane p and one point v not included in the 
plane p. In the case that a plurality of such points , ...v^ exist for the plane p, a sum A(p) of the squares of these 
points Vj (i = 1 , 2, ...m) from the plane p is expressed by Equation (5): 



.(5) 



[0053] In the data reducing method according to this embodiment, if there are n planes of a polygon model which 
are influenced by the edge shrinking, for example, when data are reduced by edge shrinking, distances of a plurality 
of original vertices (original three-dimensional shape data) involved in the data reduction processing for the approxi- 
mation from the planes are estimated for all such planes. Accordingly, the sum A(p) of the squares of the distances is 
calculated by Equation (5) for each of such planes, and e which is obtained by dividing a sum of the sums A (Pj) by a 
total number M of all vertices involved in the data reduction processing as seen in Equation (6) is assumed as an error 
estimation value: 



£ « 




[0054] Further, Equation (7) is obtained by substituting 



Kv. 'P 
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of Equation (5) for A(P/) in Equation (6): 

..^—p. ...,„ 

l"' 

10 It should be noted that in Equation (7), 

P/ denotes a column vector for an equation of a plane Pj subject to the deformation; 

Pj denotes an inverse matrix of P/ ; 

Kvy denotes an error matrix for a j-th vertex Vj involved in the deformation of the plane Pj when an error matrix for 
15 data of an original vertex v involved in the deformation of the plane Pj is defined by Kv = VV^( ^/denotes a column 

vector for the data); 

m/ denotes a sum of error matrix Kvy for a vertex Vj involved in only the deformation of the i-th plane Pj; and 
n denotes a sum of planes Pj to be deformed. 

20 [0055] In Equation (6), A(pj) is a sum of the squares of all original vertices involved in the approximation of the i-th 
plane Pj from the plane Pj, and mj is the number of terms in the sum of the squares of the distances. 
[0056] In the calculation by Equation (6), since a matrix KVj is a symmetric matrix as shown in Equation (4), only 
diagonal elements and elements above them may be stored in the memory for the matrix KV| calculated for the original 
vertex Vj. This enables a memory capacity necessary for the calculation to be reduced. Further, the respective terms 

25 rrij in the denominator of Equation (6) are the number of terms of the sum of the squares of the distances in the i-th 
approximated plane Pj. Since the element at the right bottom corner of the matrix Kv is constantly "1 " and the element 
at the right bottom corner of a matrix EKv in the calculation of the sum of the squares of the distances is mj, calculation 
by Equation (6) may be conducted using the value of this element. Thus, the calculation is easier because it is not 
necessary to separately calculate the number of terms mj. 

30 [0057] Further, it may be appreciated to use weighted error matrix for the calculation of the error e. Specifically, there 
are provided a weight factor Wj for each vertex Vj in consideration of the influence of the vertex V, to an approximated 
polygon model or other causes. The weighted error matrix is obtained by multiplying the error matrix Kv/for the vertex 

by a predetermined weight factor Wj for the vertex as w,- ■ Kvf. 
[0058] In the case that Equation (5) is applied in estimating an error for a certain plane p, if original vertices Vj involved 

35 in the processing of edge shrinking or surface shrinking for the plane p are specified, a sum EKVj of the matrices KVj 
obtained by the three-dimensional shape data of these vertices v, serves as a parameter indicating a degree of the 
error in the edge shrinking or surface shrinking. In the description below, ZKVj is referred to as an error matrix. 
[0059] Next, a procedure of the three-dimensional shape data reduction processing is described with reference to 
FIG. 7. The three-dimensional shape data of the measurement object G obtained by the three-dimensional shape 

40 measuring apparatus 1 are inputted to the shape data processor 3, in which a triangular polygon model is generated 
in the polygon generator 31 (Step #1 ). Further, an operator inputs a permissible value used to judge whether the edge 
shrinking processing can be performed via the permissible value input device 33 (Step #3). 

[0060] Subsequently, estimation values £j (j = 1 , 2, ...N, N is a total number of edges) in the case that all edges of 
the original polygon model are shrunk are calculated by Equation (6) in the estimation value calculator 321 (Step #5). 

45 For instance, in the example of FIG. 4, estimation values 8j (j = 1, 2, ...19) in the case of shrinking the edges Ed-j^ (r = 
2 to 5, 9, 10), Ed2s (s = 5, 6, ...9) and Ed^^^ (t = 3, 4, ...10, where t+1 =11 is assumed to be t+1 =3) are calculated. 
[0061] At this time, if the edge Ed-, 2 is shrunk by, for example, converging the vertices v-,, V2 to the vertex v^', the 
original planes are all deformed and eight planes p^, P2, ... p-io (hereinafter, planes after the deformation are referred 
to as approximated planes) are formed by the vertex v-,' and the vertices V3 to v-|o not involved in the edge shrinking. 

50 A sum of the squares of the distances A(P|^) (k = 1 , 2, ...8) are calculated for each of these approximated planes p^, a 
total number M of the elements included in all the sums A(Pk) is calculated, and an evaluation value e(Edi2) in the 
shrinking of the edge Ed-, 2 is calculated by Equation (6). The evaluation value e(Ed-|2) is expressed by Equation (8): 
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[0062] At this time, tine plane p-| is defined by the vertices Vj', Vg, V4. Since the vertices involved in the edge shrinking 
are only v-,, V2 if the vertices v-|, V2, V3, V4 are actually measured data, the sum A(p-|) of the squares of the distances 
is calculated by Equation (5) using an error matrix (KV1+KV2) relating to these vertices v^, V2. If, in this state, the vertex 
Vj' and another vertex v^ (unillustrated actually measured data) are converged to a vertex v^" (not shown) to approx- 

5 imate the approximated plane p^ further to an approximated plane p-,' (not shown), the vertices involved in the edge 
shrinking leading to the approximated plane p^' are v-,, V2, v^^. Thus, a sum A(Pi-) of the squares of the distances is 
calculated by Equation (5) using an error matrix (Kvi+Kv2+KVr^) relating to these vertices v^, V2, v^^. 
[0063] In the second edge shrinking, it is assumed that an error matrix Kv-,' for the vertex v-|' is (KV-1+KV2), and the 
sum A(pi') of the squares of the distances is calculated by Equation (5) using the error matrix Kv^' for the vertex v-,' 

10 and an error matrix Kv^ for the vertex v^. 

[0064] Accordingly, a general theory of the edge shrinking processing is as follows. If the vertices v^, V2, Vg, V4 are 
data generated by the edge shrinking instead of being actually measured data in FIG. 4, the error matrix Kv^' for the 
vertex v-,' defining the approximated plane p-, is Kv-|'= Kv^+Kv2 when the error matrices for the respective vertices v-,, 
V2, V3, V4 are assumed to be Kv^ , KV2, Kvg, KV4. The sum A(p^ ) of the squares of the distances is calculated by Equation 

15 (5) using an error matrix Kp-|(= Kvi'+Kv3+Kv4 = KV1+KV2+KV3+KV4) obtained by adding the error matrices Kv-|', KV3, 
KV4 for the three vertex v^', V3, V4 defining the approximated plane p-,. 
[0065] An equation used to calculate the sum A(p-|) of the squares of the distances is: 

20 MP^)-P/■Kp^•P^ (9) 

Here, P-, is a column vector of the approximated plane p^ and m-, is a total number of the error matrices Kv included 
in the error matrix Kp^. 

[0066] For the approximated planes P2, P3, ... Pg, sums A(P2), A(P3), ... A(P8) of the squares of the distances are 
25 calculated by the same method as for the aforementioned approximated plane p-,, and are substituted into Equation 

(8) to calculate the error estimation value 8(Ed-|2) in the shrinking of the edge Ed^2- 

[0067] The error estimation values e(Ed-|^), £(Ed2s), e(Ed^^^i) in the case that the other edges Ed^^ Ed2s, Ed^^^^ are 
shrunk by the same method are calculated. 

[0068] Subsequently, the error estimation values ej are sorted in the estimation value sorter 322 (Step #7) and it is 
30 discriminated whether the minimum estimation value Min{8j} is equal to or below the permissible value inputted in Step 
#3 (Step #9). If the minimum estimation value Min{£j} exceeds the permissible value (NO in Step #9), the data reduction 
processing by the edge shrinking is completed since a change of the shape resulting from application of the edge 
shrinking to all the edges is impermissibly large. On the other hand, if the minimum estimation value Min{ej} is equal 
to or below the permissible value (YES in Step #9), the edge corresponding to the minimum estimation value Min{8j} 
35 is shrunk in the edge/surface shrinking section 323 (Step #11). For example, if the minimum estimation value Min{8j} 
is 8(Ed-|2) in FIG. 4, the edge Ed-12 is shrunk as shown in FIG. 4. 

[0069] Subsequently, this routine returns to Step #7 after error estimation values (q = 1 , 2, ...) are calculated for 
the edges of the polygon model after the edge shrinking that have been influenced by the edge shrinking (Step #13). 
A minimum estimation value Min{8q} of these estimation values 8q is calculated and the edge corresponding thereto is 
40 shrunk (Steps #7 to #11). Thereafter, the edge shrinking is repeated edge by edge (loop of Steps #7 to #13), and the 
data reduction processing is completed when the minimum estimation value iVlin{eq} exceeds the permissible value 
(NO in Step #9). 

[0070] By performing the data reduction processing as above, the polygon model after the data reduction is approx- 
imated from the original polygon model such that the error estimation values from the original measurement data are 

"^5 equal to or below the permissible value. 

[0071] As described above, in the data reduction method according to this embodiment, the distances of the planes 
of the polygon model approximated by the edge shrinking or surface shrinking from the original vertices (measured 
three-dimensional shape data) are treated as the error estimation values. Thus, the data reduction processing can be 
performed by the estimation standard intuitively easily understandable. 

50 [0072] Further, since the data reduction is performed such that the evaluation values lie within the specified permis- 
sible error range, it can be performed at a high speed while relatively precisely controlling the errors of the three- 
dimensional shape data forming the polygon model after the approximation. 

[0073] Particularly, reduction in the reliability of the error evaluation values can be suppressed even in the case that 
random noise data are mixed into the measured three-dimensional shape data or data reduction is performed by 
55 converging a multitude of vertices to one point. 

[0074] FIGS. 8 and 9 show a comparison of the error estimating method according to this embodiment and the 
conventional error estimating method in the case that random noise data are mixed into the measured three-dimen- 
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sional shape data. 

[0075] FIG. 8 shows the edge shrinking of a portion (shown in section) where discontinuous surfaces are created 
on a polygon model by random noise data. If vertices v-|, V2 of the portion where large discontinuities are created on 
the outer surface of the polygon model are converged to a vertex v-,' as shown in FIG. 8, error estimation is made 
based on, for example, distances d^', 62 of the vertex v^' to a plane p-14 including vertices v^, V4 and a plane P25 
including vertices V2, V5 as shown in FIG. 9 according to the data reducing method disclosed in the above paper "Surface 
Simplification Using Quadric Error Metrics". Thus, the calculated error estimation values are smaller than sensory 
approximation errors. Contrary to this, according to the data reducing method of this embodiment, error estimation is 
based on, for example, distances d^, 62 of the vertices v^, V2 to a plane P45 including vertices v^', V4 and V5. Thus, the 
error estimation values relatively agree with the sensual approximation errors and are not calculated smaller than the 
latter. 

[0076] FIGS. 10, 11 and 12 show a comparison of the error estimating method according to this embodiment and 

the conventional error estimating method in the case that a multitude of vertices are converged to one point. 
[0077] FIG. 1 0 shows the shrinking of a zigzag surface (shown in section) defined by eight vertices v-, , V2, ... Vq of a 
polygon model. If the vertices v-|, V2, ... Vg are converged to a vertex v^' located substantially in the middle between 
the vertices Vy and Vg as shown in FIG. 1 0, error estimation is made based on distances d^', d2', ...dg' of the vertex v^' 
to the respective zigzag planes py-,, p-12, P23, P34, P45, P56' Pes original polygon model as shown in FIG. 12 

according to the data reducing method disclosed in the above paper "Surface Simplification Using Quadric Error Met- 
rics". Thus, the calculated error estimation values are larger than sensory approximation errors. Contrary to this, ac- 
cording to the data reducing method of this embodiment, error estimation is based on, for example, distances d-|, d2, ... 
dg of the respective vertex v-,, V2, ... Vq to the plane pyg including the vertices v^', Vy, Vg as shown in FIG. 11. Thus, the 
error estimation values relatively agree with the sensual approximation errors and are not calculated larger than the 
latter. 

[0078] FIG. 1 3 is a flowchart showing a modification of the procedure of the three-dimensional shape data reduction 
processing. In the procedure of the data reduction processing shown in FIG. 7, the estimation values e are recalculated 
for the edges influenced by the edge shrinking every time the edge shrinking is performed, and a next edge shrinking 
is performed based on the estimation values e. However, in the procedure of the data reduction processing shown in 
FIG. 13, the edge shrinking is applied to edges other than those influenced by the edge shrinking without calculating 
the estimation values e for the edges influenced by the edge shrinking. 

[0079] Specifically, the flowchart of FIG. 13 differs from that of FIG. 7 in that the operation in Step #13 is changed 
to Step #13a in which shrinking of the edges influenced by the edge shrinking is prohibited, and Steps #15 and #17 
are added between Step #9 and "END". 

[0080] Accordingly, common part of the flowcharts is briefly described and only different part is described in detail 
below. 

[0081] The three-dimensional shape data of the measurement object G obtained by the three-dimensional shape 
measuring apparatus 1 are inputted to the shape data processor 3, in which a triangular polygon model is generated 
in the polygon generator 31 (Step #1 ). Further, an operator inputs a permissible value used to judge whether the edge 
shrinking processing can be performed via the permissible value input device 33 (Step #3). 

[0082] Subsequently, estimation values ej (j = 1, 2, ...N, N is a total number of edges) in the case that all edges of 
the polygon model are shrunk are calculated in the estimation value calculator 321 (Step #5), and the error estimation 
values ej are sorted in the estimation value sorter 322 (Step #7). 

[0083] Subsequently, it is discriminated whether the minimum estimation value Min{8j} is equal to or below the per- 
missible value inputted in Step #3 (Step #9). If the minimum estimation value Min{ej} is equal to or below the permissible 
value (YES in Step #9), the edge corresponding to the minimum estimation value Min{8j} is shrunk in the edge/surface 
shrinking section 323 (Step #11). 

[0084] Subsequently, application of the edge shrinking processing to the edges influenced by the previous edge 

shrinking processing is prohibited (Step #13a) and then Step #7 follows. This operation is made to apply a next edge 
shrinking processing to the edges not influenced by the edge shrinking processing this time. 

[0085] Back in Step #7, estimation values eq (q = 1 , 2, ...Q where Q<N) corresponding to the edges except those for 
which the edge shrinking is prohibited are sorted, and a minimum estimation value Min{8q} is calculated. It is then 
discriminated whether the minimum estimation value Min{eq} is equal to or below the permissible value (Step #9). If 
the minimum estimation value Min{8q} is equal to or below the permissible value (YES in Step #9), the edge corre- 
sponding to the minimum estimation value Min{8q} is shrunk in the edge/surface shrinking section 323 (Step #11). 
Further, application of the edge shrinking processing to the edges influenced by the previous edge shrinking processing 
is prohibited (Step #13a) and then Step #7 follows. 

[0086] Thereafter, by repeating the above operations, a next edge shrinking processing is performed while excluding 
the edges already influenced by the previous edge shrinking processing every time the edge shrinking processing is 
performed (loop of Steps #7, #9, #11 , #13a). 
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[0087] On the other hand, if the minimum estimation value Min{8q} becomes larger than the permissible value (NO 
in Step #9) while the edge shrinking processing is repeated, it is discriminated whether the edge shrinking has been 
prohibited for these edges (Step #15). If the discrimination result is affirmative (YES in Step #15), Step #7 follows after 
error estimation values (r = 1 , 2, ...) are recalculated for these edges (Step #17). A minimum estimation value MIn 
5 {ej of the error estimation values are calculated, and the edge corresponding to the minimum estimation value Min 
{ej is shrunk (Steps #7, #9, #11). 

[0088] If the minimum estimation value l\/lin{eq} becomes larger than the permissible value and there is no more edge 
for which the edge shrinking is prohibited (NO in Steps #9, #15), the data reduction processing by the edge shrinking 
is completed. 

10 [0089] FIG. 14 shows an exemplary data reduction process by the operation procedure shown in FIG. 13. In FIG. 
14, a polygon model Mq is an original polygon model. This polygon model Mq is formed by arraying 28 vertices v^, 
V2, ...V28 in 7 rows such that the numbers of the vertices in the respective rows are (3, 4, 5, 4, 5, 4, 3) and connecting 
adjacent vertices. 

[0090] Edges el , e2, e3, e4, e5 are selected as such to be successively shrunk by the operation procedure shown 
15 In FIG. 13 to simplify the description since they are not mutually influenced in the case that the respective edges ei (i 
= 1,2, ...5) are shrunk. 

[0091] Accordingly, the data reduction process of FIG. 14 shows a case where the edges el , e2, ...e5 of the polygon 

model Mq are successively shrunk by the operation procedure of FIG. 13. Polygon models M-,, M2, M3, M4, M5 are 
approximated polygon models obtained when the polygon model Mq is successively approximated by successively 

20 shrinking the edges e1 , e2, ...e5. 

[0092] According to the procedure of the data reduction processing shown in FIG. 13, when the polygon model Mq 
is first approximated to the polygon model by shrinking the edge el by the first edge shrinking processing, edges 
formed by a converged vertex v^ and its neighboring vertices Vg, v^q, v^^, v^g, v^g, v^q, v-19, V20 (vertices located within 
adotted-line frame in the polygon model M^) are influenced. Thus, in the next edge shrinking processing, the estimation 

25 values Eq are sorted again for the edges (those located outside the dotted-line frame of the polygon model M^) excluding 
the above edges, and the polygon model is approximated to the polygon model M2 by shrinking the edge e2 cor- 
responding to the minimum estimation value Min{eq}. 

[0093] Further, since the edges defined by a converged vertex and the vertices v^, V2, V4, Vg, Vq, v^Q' ^13, Vg 
(vertices located within a dotted-line frame in the polygon model M2) are influenced in the next edge shrinking process- 
30 ing, the estimation values are sorted again for the edges (those located outside the dotted-line frame of the polygon 
model M2) excluding the above edges, and the polygon model M2 is approximated to the polygon model M3 by shrinking 
the edge e3 corresponding to the minimum estimation value Min{8s}. 

[0094] The edge shrinking is applied to the edges e4, e5 by the same method, and the original polygon model IVIq 
is approximated to the polygon model M5 and the number of data is reduced by 5 by the five edge shrinking processings. 

35 [0095] If a process for approximating the original polygon model Mq to a 6^^ polygon model Mg (not shown) is con- 
sidered, no estimation value is recalculated for the edges having been influenced by the edge shrinking according to 
the operation procedure shown in FIG. 1 3. Thus, in the example of FIG. 1 4, when the estimation value is first calculated 
for 66 edges forming the original polygon model Mq, the evaluation values are recalculated for 50 edges forming the 
polygon model M5 without calculating the estimation values again for the approximated polygon models IVI^ to M4. 

40 Therefore, the estimation values are recalculated 50 times until the 6^^^ polygon model Mg is obtained. 

[0096] On the other hand, since the estimation values are recalculating for the edges influenced by the edge shrinking 
according to the operation procedure shown in FIG. 7, the estimation values are recalculated for 16 edges influenced 
by the edge shrinking for each of the approximated polygon models M-, to M5. Thus, according to the operation pro- 
cedure shown in FIG. 7, the estimation values are recalculated 80 (16x5) times until the 6^^ polygon model Mg is 

45 obtained. Therefore, the operation procedure shown in FIG. 13 has an advantage of high-speed processing since the 
estimation values are recalculated a reduced number of times. 

[0097] FIG. 15 shows a modification of the flowchart of FIG. 13. In the flowchart shown in FIG. 15, the edges influ- 
enced by the edge shrinking are made not to be shrunk in the next edge shrinking by replacing the estimation values 
thereof by a predetermined estimation value e^^^. larger than the permissible value. 

50 [0098] Specifically, the flowchart of FIG. 1 5 differs from that of FIG. 1 3 in that Step #1 3a is changed to Step #1 3b in 
which the estimation values of the edges influenced by the edge shrinking are set at e^ax' discrimination criteria 

in Step #15 and the objects in Step #17 whose estimation values are recalculated are changed to "edges whose 
estimation values are e^ax' Steps #19 and #21 are added after Step #17a to shrink the edges whose estimation 
values are recalculated. Accordingly, only part of the flowchart of FIG. 15 different from that of FIG. 13 is described in 

55 detail below. 

[0099] When the edge corresponding to the minimum estimation value Min{8j} is shrunk in the edge/surface shrinking 
section 323 (Step #1 1 ), the estimation values of the edges influenced by this edge shrinking processing are set at the 
predetermined estimation value e^^y^ larger than the permissible value (Step #13b) and then Step #7 follows. This 
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operation is performed so as not to shrinl< tliese edges in tine next slirinking processing, i.e., to apply tine edge slirinl<ing 
only to the edges not influenced by the edge shrinking processing this time. 

[0100] Back in Step #7, evaluation values Ej (q = 1, 2, ...N-1) are sorted, and a minimum estimation value Min{ej} is 
calculated. It is then discriminated whether the minimum estimation value Mln{ej} Is equal to or below the permissible 

5 value (Step #9). If the minimum estimation value Min{8|} is equal to or below the permissible value (YES in Step #9), 
the edge corresponding to the minimum estimation value Minjej} is shrunk in the edge/surface shrinking section 323 
(Step #11). 

[0101] Subsequently, the estimation values of the edges influenced by the edge shrinking processing this time are 
set at £^3x (Step #13b) and then Step #7 follows again. Thereafter, by repeating the above operations, the next edge 

10 shrinking processing is performed while setting the estimation values of the edges influenced by the edge shrinking 
process at e^^^ every time the edge shrinking processing is performed (loop of Steps #7, #9, #11, #13). 
[0102] On the other hand, if the minimum estimation value Min{ei} becomes smaller than the permissible value (NO 
in Step #9) while the edge shrinking processing is repeated, it is discriminated whether there is any edge having the 
estimation value of e^^ax (Step #1 5a). If there are such edges (YES in Step #1 5a), the estimation values eg (s = 1 , 2, ...) 

15 are calculated for these edges (Step #1 7a). 

[0103] A minimum estimation value Min{8s} is calculated by discriminating the error estimation values 83 (Step #19). 
If the minimum estimation value Min{8s} is equal to or below than the permissible value (YES in Step #21), Step #11 
follows in which the edge corresponding to the minimum estimation value Min{8s} is shrunk. 

[0104] On the other hand, if there is no more edge having the estimation value of e^ax (NO in Step #15a) or the 
20 minimum estimation value Min{es} of the estimation values calculated for the edges having the estimation value of 
e^ax becomes larger than the permissible value (NO in Step #21 ) while the edge shrinking processing is repeated, the 
data reduction processing by the edge shrinking is completed. 

[0105] Since the shrinking of the edges influenced by the edge shrinking is prohibited in the operation procedure of 
FIG. 1 5 as in the operation procedure of FIG. 1 3 by setting the estimation values of such edges at the estimation value 
^max largei' than the permissible value, the estimation values are recalculated a reduced number of times and high- 
speed processing can be performed. 

[0106] In the procedures of the data reduction processing shown in FIGS. 13 and 15, the estimation values are not 
calculated for the edges influenced by the edge shrinking processing, i.e., these edges are not shrunk until the edge 
shrinking processing is prohibited for all edges. However, the prohibition of the edge shrinking processing may be 

30 released if a specified condition is met such as when the edge shrinking processing is repeated a specified number of 
times or when a difference between the minimum estimation value and the permissible value lies within a specified 
range, and the edge shrinking may be made possible again by recalculating the estimation values for all edges. 
[0107] Such an arrangement enables the data reduction processing to be applied to the edges for which the edge 
shrinking processing is prohibited, but the data reduction processing can be performed within a permissible error range. 

35 Therefore, data can be maximally efficiently reduced within a permissible range. 

[0108] Since the procedures of the data reduction processings shown in FIGS. 13 and 15 are designed to increase 
the processing speed by reducing the number of times the estimation values are recalculated, the estimation values 
are not limited to the errors e expressed by Equation (6). The data reduction according to the present invention is widely 
applicable even in the case that estimation values are changed by the edge shrinking or surface shrinking and need 

40 to be recalculated after shrinking, for example, those disclosed in the aforementioned paper "Surface Simplification 
Using Quadric Error Metrics" are adopted. 

[0109] Although reduction of the three-dimensional measurement data obtained by the three-dimensional shape 
measuring apparatus 1 is described in the foregoing embodiment, the present invention is not limited thereto and is 
widely applicable to three-dimensional data created by computer graphics or the like. 

45 [0110] As described above, the three-dimensional shape data are reduced by converging two or more vertices to 
one vertex to shrink the edges or surfaces. Specified estimation values are calculated for all the surfaces deformed 
by shrinking the edges or surfaces of the polygon model based on the distances between the respective surfaces and 
all the original vertices involved in the deformation of these surfaces. The three-dimensional shape data are reduced 
by edge shrinking or surface shrinking when the estimation values are equal to or below a predetermined permissible 

50 value. Accordingly, the error can be controlled by the estimation standard intuitively easily understandable, i.e., by 
approximation errors from the original three-dimensional shape data. 

[0111] Also, even in the case that, for example, the orientation of the surface locally experiences a sudden change 
into a zigzag surface in the polygon model, the data reduction processing can be suitably performed while the errors 
are relatively precisely controlled. 
55 [01 12] Particularly, since the estimation values are assumed as errors e defined by Equation (6), the errors 8 can be 
easily calculated only by storing an error matrix calculated for the original three-dimensional shape data. Thus, a high- 
speed processing is possible and a memory necessary for the calculation of the data reduction processing can be 
efficiently utilized. 
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[0113] Furthermore, the three-dimensional shape data is reduced by repeating the edge or surface based on esti- 
mation values. However, the application of the shrinking to the edges influenced by the previous edge shrinking is 
prohibited. Accordingly, recalculation of estimation values can be executed a reduced number of times, which thus 
enables a higher speed processing. 

[0114] The converged vertex portions, which are influenced by the edge or surface shrinking, are given with a pre- 
determined value to place these portions out of the shrinking, which can thus assure easier shrinking prohibition to a 
particular portion. 

[0115] Moreover, in the case that a specified condition is satisfied, the prohibition is released to allow the shrinking 
processing to the prohibited portion again. Accordingly, data reduction can be accomplished as efficiently as possible 
within a permissible range. 

[01 16] As described above, a method for reducing three-dimensional shape data, comprises the steps of: calculating 
estimation values for portions to be deformed by converging two or more vertices of a polygon model; reducing the 
number of data for the polygon model by converging two or more vertices of the polygon model when the estimation 
values are equal to or below a predetermined permissible value. 

[0117] The portions to be deformed may be surfaces of the polygon model, and estimation values for the surfaces 
are calculated based on distances between the respective surfaces and all the original vertices involved in the surface 
deformation. 

[01 18] The estimation value may be an error e defined by Equation (A): 




wherein: 



P/ denotes a column vector for an equation of a plane Pj subject to the deformation; 

P^denotes an inverse matrix of Pj ; 

Kvy denotes an error matrix for a j-th vertex Vj involved in the deformation of the plane Pj when an error matrix for 
data of an original vertex v involved in the deformation of the plane Pj is defined by Kv =V^- l/^(V^denotes a column 
vector for the data); 

/T7/ denotes a sum of error matrix Kvyfor a vertex Vj involved in only the deformation of the i-th plane Pj; and 
n denotes a sum of planes Pj to be deformed. 

[01 19] The polygon model may include a number of triangular polygons. 

[0120] The number of data for the polygon may be reduced by converging two or more vertices of one portion of the 
polygon model after another portion, repeatedly. Before each data reduction, the portion having been involved in all 
the previous data reductions may be defined as a reduction prohibition area. A succeeding data reduction may be 
applied to a portion other than the reduction prohibition area. 

[0121] The reduction prohibition area may be released if a predetermined condition is satisfied. 

[0122] The predetermined condition may be a state that there is no portion to be deformed by converging two or 

more vertices. 

[0123] The predetermined condition may be a state that a predetermined number of data reductions are completed. 

[0124] An apparatus for reducing three-dimensional shape data, comprises: estimation value calculating means for 
calculating estimation values for portions to be deformed by converging two or more vertices of a polygon model; and 
data reducing means for reducing the number of data for the polygon model by converging two or more vertices of the 
polygon model when the estimation values are equal to or below the predetermined permissible value. 
[01 25] The estimation value calculating means may calculate estimation values for surfaces to be deformed by shrink- 
ing edges or surfaces based on distances between the respective surfaces and all the original vertices involved in the 
surface deformation. 

[0126] The estimation value may be an error e defined by Equation (A): 
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...(A) 



wherein: 



P/ denotes a column vector for an equation of a plane Pj subject to the deformation; 
P/^ denotes an inverse matrix of Pj ; 

Kvy denotes an error matrix for a j-th vertex Vj involved in the deformation of the plane Pj when an error matrix for 
data of an original vertex v involved in the deformation of the plane Pj is defined by Kv=V V^ ( l/denotes a column 
15 vector for the data); 

m/ denotes a sum of error matrix Kvyfor a vertex Vj involved in only the deformation of the i-th plane Pj; and 
n denotes a sum of planes Pj to be deformed. 



[0127] The estimation value calculating means may calculate an error matrix Kv for each of all the vertices v consti- 
20 tuting an original polygon model. When two or more vertices Vp (p=1, 2, ...) are converged into one vertex Vq, the 

estimation value calculating means calculates an error e in accordance with the above-defined Equation (A) by using 

a sum of error matrix KVp for each converged vertex Vp as an error matrix KVq for the vertex Vq. 

[0128] The estimation value calculating means may include a memory which stores only diagonal elements and 

elements above them for the error matrix calculated for the original vertex Vj. The estimation value calculating means 
25 may restore the error matrix using the diagonal elements and elements above them stored in the memory, and calculate 

an error e in accordance with the above-defined Equation (A) by using the restored error matrix. 

[0129] The estimation value calculating means may calculate 




of the above-defined Equation (A) by adding an element at the right bottom corner of error matrix EKVj (j=1, 2, ... mj) 
35 for each plane Pj when the column vector V for data of the original vertex is defined by Equation (B): 



40 



. ..(B) 



45 [01 30] The estimation value calculating means may calculate an error 8 in accordance with the above-defined Equa- 
tion (A) by using a weighted error matrix Wj ■ Kvy for the vertex Vj in place of Kvj. 
[0131] The polygon model may include a number of triangular polygons. 

[0132] The data reducing means may reduce the number of data for the polygon model by converging two or more 
vertices of one portion of the polygon model after another portion, repeatedly. There may be further provided judging 
50 means for judging based on a calculated estimated value whether a portion of the polygon model is permissible for 
the deformation; and prohibition area defining means for defining a portion having been involved in all the previous 
data reductions executed by the data reducing means as a prohibition area to keep a succeeding data reduction from 
being applied. 

[0133] The prohibition area defining means may change the estimation value for an involved portion to a predeter- 
55 mined value to render the involved portion impermissible by the judging means. 

[0134] There may be further provided prohibition releasing means for releasing the prohibition area according to a 
predetermined condition. 

[0135] The predetermined condition may be a state that there is no portion to be deformed by converging two or 
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more vertices. 

[0136] The predetermined condition may be a state that a predetermined number of data reductions are completed. 
[0137] The judging means may include: minimum estimation value calculating means for calculating a minimum 
estimation value; comparing means for comparing the minimum estimation value with a permissible value; and deter- 
mining means for determining a portion having the minimum estimation value as a portion to be converged when the 
minimum estimation value is equal or below the permissible value. 
[0138] The portion to be converged may be an edge of a triangular polygon. 
[0139] The portion to be converged may be a surface of a triangular polygon. 



Claims 

1. A method for reducing three-dimensional shape data, comprising the steps of: 

calculating estimation values for portions to be deformed by converging two or more vertices (Vj) of a polygon 
model (IVIj); 

reducing the number of data for the polygon model (IVIj) by converging two or more vertices (Vj) of the polygon 
model (Mj) when the estimation values are equal to or below a predetermined permissible value. 

2. A method according to claim 1 , wherein the portions to be deformed are surfaces of the polygon model (Mj), and 
estimation values for the surfaces are calculated based on distances between the respective surfaces and all the 
original vertices (Vj) involved in the surface deformation. 

3. A method according to claim 2, wherein the estimation value is an error e defined by Equation (A): 



£ = 




wherein: 

P,- denotes a column vector for an equation of a plane p, subject to the deformation; 
P;^ denotes an inverse matrix of P-, 

Kvy denotes an error matrix for a j-th vertex Vj involved in the deformation of the plane Pj when an error matrix 
for data of an original vertex v involved in the deformation of the plane Pj is defined by Kv = V ( V^denotes 

a column vector for the data); 

/77/ denotes a sum of error matrix Kvy for a vertex Vj involved in only the deformation of the i-th plane Pj; and 
n denotes a sum of planes Pj to be deformed. 

4. A method according to any one of claim 1 to 3, wherein the polygon model (IVIj) includes a number of triangular 
polygons. 

5. A method according to claim 1 , wherein the number of data for the polygon is reduced by converging two or more 

vertices (vi) of one portion of the polygon model (Mj) after another portion, repeatedly, and before each data re- 
duction, the portion having been involved in all the previous data reductions is defined as a reduction prohibition 
area, and a succeeding data reduction is applied to a portion other than the reduction prohibition area. 

6. A method according to claim 5, wherein the reduction prohibition area is released if a predetermined condition is 
satisfied. 

7. A method according to claim 6, where the predetermined condition is a state that there is no portion to be deformed 
by converging two or more vertices (Vj). 

8. A method according to claim 6, wherein the predetermined condition is a state that a predetermined number of 
data reductions are completed. 
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9. An apparatus for reducing three-dimensional sliape data, comprising: 

estimation value calculating means (321) for calculating estimation values for portions to be deformed by 
converging two or more vertices (Vj) of a polygon model (Mj); and 

data reducing means (323) for reducing the number of data for the polygon model (M|) by converging two or 
more vertices (Vj) of the polygon model (Mj) when the estimation values are equal to or below the predetermined 
permissible value. 

10. An apparatus according to claim 9, wherein the estimation value calculating means (321) calculates estimation 
values for surfaces to be deformed by shrinking edges or surfaces based on distances between the respective 
surfaces and all the original vertices (Vj) involved in the surface deformation. 

11. An apparatus according to claim 10, wherein the estimation value is an error e defined by Equation (A): 



P/ denotes a column vector for an equation of a plane Pj subject to the deformation; 

P/^ denotes an inverse matrix of P-, 

/<vy denotes an error matrix for a j-th vertex Vj involved in the deformation of the plane Pj when an error matrix 
for data of an original vertex v involved in the deformation of the plane Pj is defined by Kv = V ( V/denotes 
a column vector for the data); 

07/ denotes a sum of error matrix Kvy for a vertex Vj involved in only the deformation of the i-th plane Pj; and 
n denotes a sum of planes Pj to be deformed. 

12. An apparatus according to any one of claims 9 to 11, wherein the estimation value calculating means (321) cal- 
culates an error matrix Kv for each of all the vertices v constituting an original polygon model (Mj); and when two 
or more vertices Vp (p=1, 2, ...) are converged into one vertex Vq, the estimation value calculating means (321) 
calculates an error e in accordance with the above-defined Equation (A) by using a sum of error matrix KVp for 
each converged vertex Vp as an error matrix KVq for the vertex Vq. 

13. An apparatus according to claim 12, wherein the estimation value calculating means (321) includes a memory 
which stores only diagonal elements and elements above them for the error matrix calculated for the original vertex 
Vj, and the estimation value calculating means (321) restores the error matrix using the diagonal elements and 
elements above them stored in the memory, and calculates an error e in accordance with the above-defined Equa- 
tion (A) by using the restored error matrix. 

14. An apparatus according to claim 11, wherein the estimation value calculating means (321) calculates 



e = 




...(A) 



n 




wherein: 



n 




of the above-defined Equation (A) by adding an element at the right bottom corner of error matrix ZKVj (j=1 , 2, 
mj) for each plane pj when the column vector V for data of the original vertex is defined by Equation (B): 
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v.. 



...(B) 



15. An apparatus according to claim 11 , wlierein tine estimation value calculating means (321 ) calculates an error e in 
accordance witli the above-defined Equation (A) by using a weighted error matrix Wj -Kvyfor the vertex Vj in place 



of Kvj. 



1 6. An apparatus according to any one of claims 9 to 1 5, wherein the polygon model (IVIj) includes a number of triangular 
polygons. 



17. An apparatus according to claim 9, wherein the data reducing means (323) reduces the number of data for the 
polygon model (Mj) by converging two or more vertices (Vj) of one portion of the polygon model (IVIj) after another 
portion, repeatedly, further comprising: 

judging means forjudging based on a calculated estimated value whether a portion of the polygon model (Mi) 
is permissible for the deformation; and 

prohibition area defining means for defining a portion having been involved in all the previous data reductions 
executed by the data reducing means as a prohibition area to keep a succeeding data reduction from being 
applied. 



18. An apparatus according to claim 17, wherein the prohibition area defining means changes the estimation value 
for an involved portion to a predetermined value to render the involved portion impermissible by the judging means. 



19. An apparatus according to claim 17 or 18, further comprising prohibition releasing means for releasing the prohi- 
bition area according to a predetermined condition. 



20. An apparatus according to claim 19, where the predetermined condition is a state that there is no portion to be 
deformed by converging two or more vertices (Vj). 

21. An apparatus according to claim 19, wherein the predetermined condition is a state that a predetermined number 
of data reductions are completed. 



22. An apparatus according to any one of claim 17 to 21 , wherein the judging means includes: 



minimum estimation value calculating means for calculating a minimum estimation value; 
comparing means for comparing the minimum estimation value with a permissible value; and 
determining means for determining a portion having the minimum estimation value as a portion to be converged 
when the minimum estimation value is equal or below the permissible value. 

23. An apparatus according to any one of claims 1 7 to 22, wherein the portion to be converged is an edge of a triangular 
polygon. 



24. 



An apparatus according to any one of claims 17 to 23, wherein the portion to be converged is a surface of a 
triangular polygon. 
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FIG.3 
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FIG.4 
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FIG.7 
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FIG.13 
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FIG.15 
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